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ABSTRACT 

We use photometric and spectroscopic observations of the eclipsing binaries 
V65, V66 and V69 in the field of the globular cluster M4 to derive masses, radii, 
and luminosities of their components. The orbital periods of these systems are 
2.29, 8.11 and 48.19 d, respectively. The measured masses of the primary and 
secondary components (M p and M s ) are 0.8035±0.0086 and 0.6050±0.0044 M 
for V65, 0.7842±0.0045 and 0.7443±0.0042 M for V66, and 0.7665±0.0053 and 
0.7278±0/0048 M for V69. The measured radii (R p and R s ) are 1.147±0.010 
and 0.6110±0.0092 R Q for V66, 0.9347±0.0048 and 0.8298±0.0053 R Q for V66, 
and 0.8655±0.0097 and 0.8074±0.0080 R for V69. The orbits of V65 and V66 
are circular, whereas that of V69 has an eccentricity of 0.38. Based on systemic 
velocities and relative proper motions, we show that all the three systems are 
members of the cluster. We find that the distance to M4 is 1.82±0.04 kpc - in 
good agreement with recent estimates based on entirely different methods. We 
compare the absolute parameters of V66 and V69 with two sets of theoretical 
isochrones in mass-radius and mass-luminosity diagrams, and for an assumed 
[Fe/H] = -1.20, [cu/Fe] = 0.4, and Y = 0.25 we find the most probable age of M4 
to be between 11.2 and 11.3 Gyr. CMD- fitting with the same parameters yields 
an age close to, or slightly in excess of, 12 Gyr. However, considering the sources 
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of uncertainty involved in CMD fitting, these two methods of age determination 
are not discrepant. Age and distance determinations can be further improved 
when infrared eclipse photometry is obtained. 

Key words: binaries: close binaries: spectroscopic globular clusters: indi- 
vidual (M4) stars: individual (V65-M4, V66-M4, V69-M4) 
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individual (M4) - stars: individual (V65 M4, V66 M4, V69 M4) 



1. Introduction 

Detached eclipsing double-line binaries (DEBs) are the primary source of the observa- 
tional data concerning stellar masses and radii. When supplemented by luminosities derived 
from parallaxes, empirical relations between color and effective temperature, or fits to dis- 
entangled spectra, they enable fundamental tests of stellar evolution models. For many 
field Population I binaries with components at solar mass or larger, modern high accuracy 
measurements of masses, r adii and luminosities are in general agreemen t with theoretical pre- 
dictions (see, for example, lLacy et al.ll2005l 120081 ; IClausen et al.ll2008l). Similar encou raging 
result s are obtained for bin aries in the old o pen clusters NGC 188 (IMeibom et al.ll2009r). NGC 



2243 (IKaluznv et alJl2006h . and NGC 6791 (IGrundahl et alJl2008 



Brogaard et al.ll201ll ). On 



the other hand, the models seem to underestimate the radii of numerous K and M dwarfs in 



in 



Blake et al. 


(2008). 


Torres et al. 


(2010 


), and 


Kraus et al. 


(2011) 



The situation is less clear for Popu lation II stars, for whic h only a very few DEBs with 
main-sequence components are known ([Thompson et al.l 120101 ). and there is an urgent need 
to locate and study such systems. Within the series Cluster AgeS Experiment (CASE), this 
is the second paper devoted to the study of globular-cluster DEBs with main-sequence or 
subgiant components. The general goal of CASE is to determine the basic stellar parameters 
(masses, luminosities, and radii) of the components of cluster binaries to a precision better 
than 1% in order to m easure cluster ages and distances, and to test stellar evolution models 
(IKaluznv et al.l 120051 ) . The meth ods and assumpti ons we employ utilize basic and simple 
approaches, following the ideas of iPaczyhskil (119971 ) and lThompson et al.l (120011) . Previous 
CASE papers analyzed blue straggler systems in uj Cen ( Kaluzny et al.l2007a ) and 47 Tuc 



*This paper includes data gathered with the 6.5-m Magellan Baade and Clay Telescopes, and the 2-5 m 
du Pont Telescope located at Las Campanas Observatory, Chile. 
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kaluznv etaDl2007bh . an SB1 binary in NGC 6397 faaluznv et "al] l2008h and the binary 



V69- 47 Tuc, which is an SB2 system with main-sequence components ([Thompson et al. 



20101 ) . 



The present paper is devoted to the analysis of three DEBs, V65, V66 and V69, all 
members of the globular cluster M4. We use radial velocity and photometric observations 
to determine accurate masses, luminosities, and radii of the components of these systems. 

Section H] describes the photometric observations and the determination of orbital ephe- 
merides. Section [3] presents the spectroscopic observations and the radial- velocity measure- 
ments. The combined photometric and spectroscopic solutions for orbital elements and 
component parameters are obtained in Section HJ while the cluster membership of the three 
DEBs is discussed in Section [5J In Section [5J we compare the derived parameters to a se- 
lection of stellar evolution models, with an emphasis on estimating the age of the system. 
Finally, in Section [7] we summarize our findings. 



2. Photometric observations 

Our survey for eclipsing binaries in M4 began in July 1995 with a 2-week observational 
campaign at CTIO. We monitored the cluster in the B- and V-bands, using the 2K 2 TEK2 
camera attached to the 0.9-m telescope. A single eclipse of V66 was detected, occurring 
on 1995 July 18 UT. The survey continued on the 1.0-m Swope telescope at Las Campanas 
Observatory (LCO) using five different CCD cameras (Ford, TEK3, SITel, SITe2, and SITe3) 
over several seasons during the period 1996 - 2009. Most of the observations were obtained 
with the SITe3 camera, with a plate scale of 0.435 arcsec/pixel and a field size of 14. 8x22.8 



arcmin . Some early results from this survey were presented in iKaluzny et al.l (119971 ). 



The first eclipses of V65 and V69 were detected on 1996 April 16 UT and 1996 April 
21 UT, respectively. In the same year preliminary but reliable ephemerides for V65 and V66 
were determined. A part of a secondary eclipse of V69 was observed on 1998 August 17 
UT. However, because of the relatively long period for this system, the data were still insuf- 
ficient to establish a unique ephemeris. An initial ephemeris was derived from constraints 
provided by substantial out-of-eclipse photometry performed in 1996 - 2005, and radial ve- 
locity observations obtained with the MIKE spectrograph on the Magellan Baade and Clay 
telescopes (see Section [3]). From 1998 until 2009 we also observed our targets with the 2.5-m 
du Pont telescope at LCO equipped with the TEK5 2K 2 camera with a plate scale of 0.259 
arcsec/pixel. Between 2001 and 2009, several eclipses were covered for V65 and V69, and 
some out-of-eclipse data were collected for V69. In all observing runs, on each telescope, the 
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same B and V niters were used. Times of minima are presented in Tables HJ [21 and [3j 



Linear ephemerides provide adequate fits to the photometric data for all three binaries, 
these are given by Equations (1). In the case of V65 and V66 the ephemerides were obtained 
from the moments of minima calculated from individual light curves. This procedure could 
not be applied to V69, as in only three eclipses both the descending and the ascending branch 
was observed. In the cas e of V69, the period was derived using the algorithm developed by 
Lafler fc Kinmannl (119651 ) , and the moment of the primary eclipse was found from the phased 
l ight curve with an improved version of the method developed by iKwee fc van Woerden 
( 119561 ). For V65, only du Pont observations were used, reduced with the image subtraction 
technique. For V66, we also made use of photometry measured with profile fitting in images 
collected with the Swope and 0.9-m CTIO telescopes. 



tt j pi V65 



HJD 



V69 



2449905.49577(43) + 2.29304564(26) x E 
2449900.41965(26) + 8.11130346(85) x E 
2450048.34890(14) + 48.1882687(6) x E 



Table [4] gives the equatorial coor dinates of the three y ariables analysed in this paper. 
They are tied to the UCAC3 system (jZacharias et al. 2010 ) for V65 and V66 , and to the 



GSC-1.0 system (e.g. lLasker et al.l Il990l ) for V69. Finding charts prepared from du Pont 



TEK5 \^-band images are shown in Figure [TJ 

V65 is blended with two stars located at angular distances of 0.37" and 0.57". The 
variable is the brightest component of the blend . All three stars are in cluded in ACS-HST 
HST F606W/F814W photometry published bv knderson et al.l koO$ \. with star identifi- 
cation numbers 13362, 13372 (V65), and 13374. The HST photometry is listed in Table 
We determined the rectangular coordinates of V65 and its two close visual compan- 
ions on the refe r ence i mage by transforming their positions from ACS/HST photometry of 



Anderson et al.l (120081 ) . and the presence of the comp anions was taken into account while 



extracting light curves of this binary. V66 is listed by I Anderson et al.l (120081 ) as star 5795. 
It does not suffer from any blending which could affect the ground based photometry. V69 
is not present in any of the currently available HST images of M4. The object is located at 
an an gular distanc e of 7.69' from the center of the cluster, beyond the half-light radius of 
4.33' (jHarrid Il996l . 2010 edition). The du Pont images show no evidence of any unresolved 
visual companions to V69. 
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2.1. Light curves and calibration of photometry. 



The B and V light curves of V65 and V66 were determined entirely from the du 
Pont data. They were extracted with the image subtraction technique, using methods 
and codes described in iKaluzny et al.l (I2O1O0 . In short, w e used a combi nation of ISIS 
2.1 jAlard fc Luptonl Il998h . Allstar/Daophot JstetsorJ [l987h and Daogrow JstetsorJ [l99oh 
codes, supplemented with some iraf| tasks. For these two variables, we also extracted 
light curves using profile fitting software to make sure that the image subtraction technique 
is free from systematic errors (in particular, we checked that the curves resulting from pro- 
file photometry and image subtraction have the same amplitudes). As expected, the scatter 
produced by the image subtraction technique was significantly smaller. 

B and V light curves of V69 were obtained from the data collected with the Swope 
telescope and the SITe3 camera. The light curves were measured using profile fitting software 
since this variable is located in a sparsely populated field and image subtraction offers no 
advantage in comparison with the traditional profile photometry. We also collected a few 
frames containing V69 with the du Pont telescope and used th ese to me a sure t he magnitude 
and color of this variable at maximum light. Observations lLandoltl ( 119921 ) photometric 
standards obtained with the du Pont telescope enabled us to transform the light curves of 
all three variables to the standard BV system. 

Full details of the procedure employed to perform the photometric measurements will be 
reported in a separate paper devoted to the photometry of other variables in M4 (Kaluzny et 
al., in preparation), along with the measurement of differential and global reddening in the 
field of the cluster. Table lists magnitudes and colors of V65, V66 and V69 at maximum 
light, together with the reddening towards each of the variables. The errors include internal 
and external uncertainties, and the last column gives the total reddening as determined by 
Kaluzny et al. (in preparation; the reddening was found by comparing turnoff colors of M4 
and the low-extinction cluster NGC 6362). 

The light curve of V65 is unstable (see Figure [2]). The fast rotation of the components, 
implied by the relatively short orbital period of 2.29 d, apparently induces a strong magnetic 
activity in at least one of the stars. A clear sign of such an activity is the X-ray emission: 
the system was listed by iBassa et al.l (120041 ) as the X-ray source CX 30 with L x = 2.6 x 



10 29 erg s 1 . For the present analysis, we selected observations collected between 2008 June 
07 and 2009 June 30. During that period the light curve was flat between the eclipses, 



X IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
AURA, Inc., under cooperative agreement with the NSF. 
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the eclipses were symmetric (note the totality of the secondary eclipse), and the binary 
was brighter than in the other observing seasons. All these observations indicate that the 
magnetic activity of the system was significantly lower than during the other observing 
seasons. 

The remaining two light curves are stable. While the eclipses of the star V66 are 
symmetric and separated by a half of the orbital period, the orbit of V69 is clearly eccentric, 
and the secondary eclipse occurs at phase 0.609. The final light curves adopted for the 
analysis of the three systems are shown in Figure [31 



3. Spectroscopic observations 



The spectra were taken with the MIKE echelle spectrograph (IBernstein et al.l 120031 ) on 
the Magellan Baade and Clay 6.5-m telescopes, using a 0.7" slit, which provided a resolution 
R ~40,000. A typical observation consisted of two 1800-second exposures of the target, 
flanking an exposure of a thorium-argon hollow-cathode lamp. A few of the exposures were 
shorter, depending on the observing conditions. The raw spectra were re duced w i th th e 
pipeline software written by Dan Kelson, following the approach outlined in iKelsonl ( 120031 ). 
The IRAF package ECHELLE was used for the post-extraction pr ocessing. The ve l ocitie s 
were measured using software based on the TODCOR algorithm (IZucker fc Mazehl Il994l ). 
kindly made available by Guillermo To rres. For velocity tem plates, we used synthetic echelle- 
resolution spectra from the library of ICoelho et al.l (120061 ) . These were interpolated to the 
values of log g and T e $f derived from the photometric solu tion (see Section HI) and assuming 
[Fe/H ] = —1.2 with an a-element enhancement of 0.4 (jCaretta et al.l 120091 ; iDotter et al. 
20101 ). The results of the velocity measurements are insensitive to minor changes in these 
parameters. The templates were Gaussian-smoothed to match the resolution of the observed 
spectra. In the case of V65, a rotational broadening was additionally applied. 

For both V66 and V69, the cross-correlations with template spectra were performed 
independently on wavelength intervals 4120-4320 A and 4350-4600 A, covering the region 
of the MIKE blue spectra with the best signal-to-noise ratio while avoiding the H 7 line. 
The final velocities adopted for the analysis were obtained by averaging the results of these 
two measurements. In the case of V65, the spectrum was contaminated by one of the stars 
blended with the target (ID 13 374 in Table [71), so that a three -dimensional extension of the 
original TODCOR algorithm ( IZucker. Torres fc Mazehl Il995l ) had to be used. The cross- 
correlation was performed on the interval 4000-4840 A to maximize the signal from the very 
faint secondary component. The second contaminating star (ID 13362 in Table |7J), which is 
fainter and more distant on the sky from the target, was not detected in the velocity cross 
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correlation functions. 

The results of velocity measurements are presented in Tables [HI and which list 
heliocentric Julian dates at mid-exposure, velocities of the primary and secondary compo- 
nents, and orbital phases of the observations, calculated according to the ephemerides given 
by Equations (1). In the case of V65 the velocities of the contaminating star are also given. 
The observed velocity curves were fit with a nonlinear least squares solution, using code 
kindly made available by Guillermo Torres. Both the observed curves and the fitted ones are 
shown in Figure HJ and the derived orbital parameters are listed in Table [TD1 together with 
errors as returned by the fitting routine. Table [TU1 also lists the velocity standard deviations 
a p and a s of the orbital solution which are a measure of the precision of a single velocity 
measurement. In all cases, the fits adopted the periods and times of primary eclipse as given 
in Equations (1). 



4. Light curve analysis and system parameters 



The analysis of the light curves was performed with the PHOEBE implementation 
flPrsa &: Zwitterll2005l ) of the Wilson-Devinney (WD) model (IWilson fc Devinneylll97ll ; IWilson 
1979). The PHOEBE/WD package utilizes the Roche geometry to approximate the shapes 



of the stars, uses Kurucz model atmospheres, treats reflection effects in detail, and, most 
importantly, allows for the simultaneous analysis of B and V data. The resulting geometri- 
cal parameters are largely determined by the higher- quality V data, while the B data serves 
mainly to estimate of the luminosity ratio L s /L p in that band. To find the best initial 
parameters for PHOEB E iterations, we solved for the V data using the JKTEBOP code 
( iSouthworth et al.ll2004l and references therein), which, unlike PHOEBE, can deal with a 
single light curve only, but it is capable of a robust search for the global minimum in the 
parameter space. 

Before solving for the light curves, it was necessary to estimate the effective temperature 
of the primary, T p . T o that end, we used the (B — V) and E{B — V) values from Table El 
and the calibration of ICasagrande et al.l ( 120101 ) . Apart from T p , PHOEBE needs to be given 
a metallicity, albedo, and gravity darkening coefficient. The user also has to specify which 
limb darkening law is to be used. We adopted a metallicity of [Fe/H] = —1.20 (see Section 
Eland Section E]). Bolometric albedo and gravity darkening coefficients were set to values 
appropriate for stars with convective envelopes: A = 0.5, and g = 0.32 (we note that all 
three systems are well detached, so that effects of reflection and gravity only weakly affect 
their light curves). We used a linear approximation for the limb darkening. Theoretical 
limb darkening coefficients in the BV bands were interpolated from the tables compiled by 
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Claret ( 2000 ). using the jktld code. 2 Q The parameters of each binary were found iteratively 



according to the following procedure: 

1. Solve for the velocity curve, as explained in Section [3) Find a preliminary solution of 
the flight curve using JKTEBOP. Feed the obtained parameters into PHOEBE. 

2. Solve for the light curves, fitting orbital inclination i, effective temperature of the sec- 
ondary T s , gravitational potentials at the surface of the primary Q p and the secondary 
Q s , and relative luminosities Lf /Lp 3 and L v s jl)^ in B and V bands. 

3. Based on the relative luminosities obtained in Step 3, calculate th e (B — V) n index 
of the primary, and update T p , using, as before, the calibrations of ICasagrande et al. 



(l2010h 



4. Repeat steps 2 and 3 until iterations converge (i.e. until the last iterated corrections 
to the parameters became smaller than the formal errors of those parameters). 

The eccentricity and argument of periastron for V69 were found from the velocity curve. We 
obtained a preliminary photometric solution in which they were allowed to vary, but since 
they changed by only 10% of the errors given in Table [TUl we decided to keep them constant 
during proper iterations. 

The residuals of the fits are shown in Figure [5] and the final values of the iterated 
parameters are given in Table [TTJ We checked that the luminosity ratios and relative radii 
of the components were practically insensitive to changes in effective temperature of the 
primary: (L s /L p )b, (L s /L p )v, r s , and r p , all changed by less than 0.3% for a ±150 K change 
in T p . The standard errors of the six parameters iterated upon by PHOEBE were found using 
a Monte Carlo procedure written in the P HOEBE-scripter, and sim ilar to that outlined in 



the description of the JKTEBOP code (see ISouthworth et al.l 12004 and references therein). 



Briefly, the procedure replaces the observed light curves B Q and V with the fitted ones Bf 
and Vf, generates Gaussian perturbations 5Bf and SVf such that the standard deviation of 
the perturbation is equal to the standard deviation of the corresponding residuals shown in 
Figure and performs PHOEBE iterations on Bf + 8Bf and Vf + SVf. Each Monte Carlo 
run produced 15000 points. Examples of the Monte Carlo diagrams are shown in Figure [6j 
The combination of spectroscopic and photometric solutions yielded absolute parameters 
of the three DEBs which we list in Table [121 The measured masses of the primary and 
secondary components (M p and M s ) are 0.8035±0.0086 and 0.6050±0.0044 M for V65, 



2 The code is available at 



http://www.astro.keele.ac.uk/jkt/codes/jktld.html 



-9 - 



0.7842±0.0045 and 0.7443±0.0042 M for V66, and 0.7665±0053 and 0.7278±0/0048 M 
for V69. The measured radii (R p - R s ) are 1.1470±0.0104 and 0.6110±0.0092 R for V66, 
0.9347±0.0048 and 0.8298±0.0053 R for V66, and 0.8655±0.0097 and 0.8074±0.0080 R for 
V69. The measured accuracy of the mass determinations ranges from 0.6% for the primary 
of V66 to 1.1% for the secondary of V65. The measured accuracy of the radii determinations 
ranges from 0.5% for the primary of V66 to 1.5% for the secondary of V65. The location 
of the components on the CMD of M4 is shown in Figure [7J This CMD shows data for 
a 3.5x3.5 arcmin 2 field whose center is located 1.4 arcmin from the center of the cluster. 
This field was chosen because it shows relatively uniform extinction. The photometry of all 
stars in Figure [7J (including the three DEBs) has been corrected for differential reddening 
(Kaluzny et al., in prep.). No attempt was made to remove nonmembers of M4. 



5. Membership of the target binaries and the distance to M4 



Upon averaging the subtracted images on a seasonal basis, we found no evidence of 
bipolar residuals at the positions of our targets. Such residuals are observe d for objects with 
notic eable proper motions with respect to the surrounding stellar field (lEyer &: Wozniak 
200ll ). Given the time base and the pixel scale of our images, we rule out motions in excess 
of 10 mas/y. Since on the proper-motion diagram of M 4 most of the field stars are separated 
from the cluster population by more than 15 mas/y ( jZloczewskil 120121 ). this is consistent 
with all the three targets being members of M4. 



The radial velocity of M4 is equal to 70.29±0.07 km s _1 ( iSommariva et al.ll2009l ). while 
the ve locity dispersion is 3 .5±0.3 km s _1 at the core, dropping marginally towards the out- 
skirts ( iPeterson et al.lll995l ). Thus, V65 and V69 are unquestionable radial- velocity members 
of the cluster. The velocity of V66 is ~2cr larger than the cluster mean, but the difference is 
too small to exclude membership. An additional argument in favor of cluster membership is 
the location of all of the components of the three DEBs on or very close to the main sequence 
of the cluster (see Figure [7J . 

Other authors have derived the distance to M4 usi ng the Baade-Wesselink method 
( jLiu fc Janesl Il990h . by astr ometrv (IPeterson et al.lll995l ). or by fitting the subdwarfs to 
the cluster's main sequence (IRicher et al.lll997l ). Their res ults are remarkab ly consistent: 
1.72±0.01, 1.72±0.14 and 1.73±0.09 kpc, respectively (see IRicher et al.l l2004h : however to 
achieve this consistency it was necessary to replace the standard yalue of the ratio R y = 
A V /E(B - V) = 3.1 with a significantly larger one (R v ~ 3.8: IRicher et al.l 20ol). A 
thorough recent study of the reddening law in the field of M4 (IHendricks et al.ll2012l ) indicates 
that, due to the intervening Scorpius-Ophiuchus dark clouds, the appropriate value is Ry = 
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3.76 ± 0.07. Their estimate of the distance to the cluster, obtained from fitting of the ZAHB 
l/-band magnitude to models, is 1.80±0.05 kpc. 



Following iHendricks et al.l ( 120121 ). we adopted Ry = 3.76, and corrected the observed 
magnitudes for extinction using the E(B — V) values from Table [6j To obtain absolute mag- 
nitudes in the V^-band we used bolometric luminosities from Table [121 ^v,& = 4.81 ± 0.03 



( [Torres! 120101 1 . and theoretical bolometric corrections obtained from models described in 
Section [6J Having the corrected observed magnitudes and absolute magnitudes, we derived 
distance moduli separately for each component of our three DEBs. The moduli are listed in 
Table [12] together with their errors which arise fro m uncertainties o f app arent magnitudes 
from Table [TTT reddening from Table [61 Ry from IHendricks et al. (l2012f ). and Mv,q from 
Torres! (120 lOh . For reasons discussed in Section [61 the effective temperatures of V65 com- 
ponents may be biased, resulting in luminosity errors that are hard to account for. Thus, 
in principle, the distance calculated from V66 and V69 only should be more reliable than 
with V65 included. With this restriction, we obtain a weighted mean distance (the weights 
being equal to in yerse errors squa r ed) o f 1.82±0.04 kpc - in excellent agreement with the 
recent estimate of IHendricks et al.l ( 120121 ). We note that when R v = 3.1 is used, the distance 
increases to 1.95 kpc, in disagreement with other estimates. Thus, our results provide an 
independent confirmation of the atypical reddening law in the field of M4. Using all six 
moduli from Table [121 one gets 1.85±0.03 kpc - a value compatible with that derived from 
V66 and V69 alone. 



6. Isochrone age analysis 



The fundamental parameters derived from the binary systems reported in this paper 
allow us to derive the ages of the individual stars, as well as an aggregate age for the 
cluster, using theoretical isochrones. However, before the model comparison is performed 
it is necessary to review the information available on the chemical composition of M4 since 
the model-based ages are sensitive to the adopted values of helium abundance, [Fe/H] and 
[a/Fe]. 

The proxi mity of M4 h a s mad e it a frequent subject of spectroscopic investigations. As 
summarized by llvans et al.l ( 119991 ). [Fe/H] determinations based on high resolution abun- 
dance analyzes up until that time range from —1.3 to —1.0 with a mean of —1.15. The 
following is a brief review of the results presented in large-scale spectroscopic surveys of M4: 



Based on the analysis of 23 stars with high resolution spectra, llvans et al.l (1l999h 
derived a mean [Fe/H] = — 1.18 ± 0.02. Their measurements of a-capture elements 
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imply a mean [a/Fe]=+0.35 (mean of 0, Mg, Si, Ca, and Ti; note also that and Mg 
exhibit significant star-to-star variations). 



Marino et all (12008( 1 derived [Fe/H] = -1.07 ± 0.01 and [a/Fe]=+0.39 ± 0.05 from high 



resolution spectra of 105 stars. 



Caretta et al.l (120091 ) derived [Fe/H] = -1.200 ± 0.053 from 103 stars observed at the 



VLT with GIRAFFE (or -1.168 ± 0.066 from 14 stars observed with UVES) where 
the quoted error is the quadrature sum of statistical and syste matic errors. Furthe r 
measurement of these same spectra yields a mean [a/Fe]=+0.51 (jCarretta et al.ll2010l ). 



Villanova fc Geislerl ( 1201 ll ) measured [Fe/H] = —1.14 ± 0.02 for 23 red giant branch 



stars located below the RGB-bump. 



For six blue horizontal branch stars. I Villanova et al.l (120121 ) determined [Fe/H] = — 1.06± 
0.02. This value holds for the supposed subpopulation of He-enriched stars with 
Y=0.29±0.01. 



We adopt a fiducial composition of [Fe/H] = — 1.2, [a/Fe]=+0.4 and Y=0.25, but will consider 
a range of these parameters while deriving the ages of the binary system components. For the 
following age analysis, we use two sets of theoretical isochrone s that are represent ative of the 
state of the art for low-mass, metal-poor stars: Dart mouth (IDotter et al.l 120081 . henceforth 
DSED) and Victoria- Regina (IVandenBerg et al.ll2012l . henceforth VR). A comparison of the 
derived stellar parameters with 10, 11, 12, and 13 Gyr DSED and VR isochrones obtained 
for the fiducial composition is shown in Figure [8] in the (M-R) and (M-L) planes. While the 
two sets of models are qualitatively similar, the quantitative differences between them lead 
to differences in derived ages. However, such differences are smaller than the uncertainties 
imposed by the observational errors. 

It is immediately clear from Figure|8]that the lowest-mass star, i.e. the secondary of V65, 
is larger in radius than either set of isochrones predicts, but that its luminosity is consistent 
with an (essentially) unevolved main sequence star. This finding is in qualitative agreement 
with the results discussed in Section [T] for nearby field binaries with similar masses, and is 
not unexpected given the dynamical and X-ray properties of V65 (see Section |2TT|) . The next 
point to notice is that the primaries of V66 and V69 yield age estimates that are consistent 
with each other in both planes. The secondaries in these systems appear to favor older ages, 
though not at a statistically significant level, particularly in the (M-L) plane. This may be 
a consequence of the anticorrelation of the radii in V66 and V69 illustrated in Figure [6J an 
overestimate of the secondary's radius implies an underestimate of the primary's one. If that 
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is the case, then the age discrepancy should more-or-less cancel out when the average age of 
each system is considered. 

In order to formally incorporate the observational uncertainties into the age analysis, 
we evaluate the age of each star on a dense grid of points within that star's 3-cr error box 3 1 in 
both the mass-radius and mass-luminosity planes. Each point at which the age is determined 
has a weight w — (1 + <5) -1 where 5 represents the distance of a point from the best value in 
(M-R) or (M-L) plane in units of the standard deviation derived from the observations. Thus 
defined, < w < 1. The ages and weights are used to construct weighted age histograms for 
each star. 

The resulting histograms for the components of V66 and V69 are displayed in Figure [9] 
V65 is omitted from the figure because its properties make it unsuitable for comparison 
with standard stellar evolution models. For completeness, the mean and standard deviations 
derived from the distributions are summarized for all stars in Table [13j As already remarked, 
the DSED and VR isochrones give ages that agree to within one standard deviation in every 
case. 

Further age uncertainties are caused by sensitivity to chemical com position and inherent 



uncer tainties in the stellar evolution model physics (a ~3% effect, see IChaboyer fc Krauss 



20021 ). We analyze the sensitivity to chemical composition in the (M-R) plane only because 
these quantities do not depend on the adopted composition, whereas the luminosity depends 
on the composition via the effective temperature. We calculate age differences (A-age) 
with respect to the age derived assuming the fiducial composition. A positive A-age value 
indicates that the model with varied composition yields an older age than the model with the 
fiducial composition. The numbers presented in Table [H] are averaged over the components 
of V66 and V69; there is a slight sensitivity of A-age to stellar mass but it is less than 
0.1 Gyr. Increasing Y decreases the age derived from the (M-R) plane while increasing 
[Fe/H] increases the age. Increasing the [a/Fe] ratio also increases the derived age, but 
the amount varies because of the way that a-enhancement is defined in each set of models. 
The Dartmouth models use a constant enhancement of the a-capture elements whereas 
the Victoria-Regina models employ an observational ly-motivated enhancement (the £ GSC 



heavy element mixture, see IVandenBerg et al.l 120121 ). The age difference between scaled- 



solar ([a/Fe]=0) and the a-enhanced mixture dep ends in detail on the am ount to which 



certain elements (most notably O) are enhanced, see I VandenBerg et al.l ( 120121 ) for a thorough 
discussion. 



3 A 3-cr error box represents the best compromise between fully sampling the (assumed normal) age 
distribution and remaining within the parameter space covered by the isochrone grids. 
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Potter et al.l ( 120091 ) discussed the insights that may be gained by comp aring stellar ages 



derived from fitting the mass-radius relation of the binary V69 in 47 Tuc ( [Thompson et al. 



20101 ) with those derived from fitting isochrones to the cluster CMD. In particular, those 
authors showed that the while the mass-radius diagram is sensitive to all aspects of the 
composition considered above, the CMD is largely insensitive to variations in Y. Furthermore, 
they found that the mass-radius and CMD ages respond differently to changes in [Fe/H]: 
while age and [Fe/H] are correlated in the (M-R) or (M-L) diagram, they are anticorrelated 
in the CMD. 

It is therefore of some value to consider the implications of the age analysis presented 
in this section to the comparison of isochrones with the CMD of M4. Figure [10] plots the 
Dartmouth models with the fiducial composition and ages of 10, 11, 12, and 13 Gyr (the 
same as shown in Figure [5}. To adjust the isochrones to the observed CMD, we adopted a 
true distance modulus of 11.34 (this is the weighted mean from six moduli listed in Table [T2|) . 
A V =1A7 and E ( B - V )=0.39. A v is the product of E(B - V) and R v = 3.76 taken from 



Hendricks et al.l ( 120121 ) . while E(B — V) itself is taken from Kaluzny et al. (in preparation), 
who derived it for a reference region with uniform reddening, in which all stars shown in 
FigureslTlandnUlreside. This derivation is based on two independent sets of observations from 
2002 and 2003. For each season they selected the best photometric night, during which over 
50 measurements of Land olt standards w ere made. The agreement with ^/-band magnitudes 



of M4 stars published by IStetsonl ( 120001 . 2012 CADC online edition) was excellent, however 
the measured (B — V) was on the average 0.023 mag larger than that of Stetson (0.018 mag 
and 0.028 mag, respectively, for the 2002 an d 2003 seasons), causing an analogous increase 



in the derived E(B — V). iHendricks et al.l ( 120121 ) obtained a slightly lower reddening of 



0.37 mag. While small differences in the photometric calibration or in the areas selected 
for the analysis can easily account for this discrepancy, the latter value is inconsistent with 
the overall agreement shown in Figure [TUl The isochrone comparison shown in Figure [TU] is 
consistent with an age > 12 Gyr, higher than derived from the binaries. We further discuss 
the age derived for M4 in Section [71 

We define the aggregate age of M4 as an average of the ages of four stars forming 
V66 and V69 systems. We calculated two averages. The first is a standard weighted one, 
including data from fits in both the (M-R) and (M-L) planes. It is equal to 11.10±0.26 and 
11.23±0.27 Gyr, respectively, for DSED and VR isochrones. The second average is obtained 
from (M-R) fits only, and in two steps. In the first step the mean age of each system is 
found by averaging the ages of the components (as explained above, this procedure removes 
effects resulting from the anticorrelation of the stellar radii). In the second step, a weighted 
average of the ages of the two systems is calculated, with the weights equal to inverse errors 
squared. For DSED and VR isochrones this procedure yields 11.25±0.42 and 11.30±0.44 
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Gyr, respectively, i.e. values entirely compatible with those obtained with the first method. 
The final (and conservative) age estimate is given by the largest range of ages resulting from 
the both methods: we may say that M4 is older than 10.8 Gyr, but younger than 11.7 Gyr, 
with the most probable age between 11.2 and 11.3 Gyr. 



Following iThompson et all ( 120101 ). and accounting for the age sensitivities listed in Table 
14, we adopt a systematic error of 0.85 Gyr arising from a 0.1 dex uncertainty in each 
of [Fe/H] and [alpha/Fe]. Based on DSED isochrones, our formal age estimate for M4 
derived from the study of the binary stars V66 and V69 is 11.25 ±0.42 ±0.85 Gyr. We note 
parenthetically that the He abundance of M4 of Y = 0.29 +/-0.01 measured by Villanova 
et al. (2012) in six blue HB stars implies an age of approximately 8 Gyr (see Table 14). A 
second burst of star formation occurring after a such a long delay seems very unlikely. 



7. Discussion and summary 



We have derived absolute parameters of the components of V65, V66 and V69 - three 
detached eclipsing binaries located on the main sequence of the globular cluster M4. The 
accuracy of our mass and radii measurements is better than 1.5% for V65, and better than 
1% for the remaining two DEBs. The reason for the lower accuracy of the parameters of V65 
is the high activity of this system which causes its light curve to be strongly variable (see 
Figure [2]). V65 is a fast rotating, X-ray active binary whose components are most probably 
puffed-up due to the presence of starspots and/or magnet ic fields, as it is often observed 
in short-period binaries of spectral type K and M (see e.g. iTorresI 120101 ). These properties, 
while interesting by themselves, make this object unsuitable for analyzes based on isochrone 
fitting. We note that chromospherically active components of eclipsing binaries appear to 
be larger a nd cooler than inact ive single stars of the same mass, but they have a similar 
luminosity ( jMorales et all 120081 ). This naturally explains why the primary of V65 is (and 
the secondary may be) located to the red of the main sequence of the cluster. 

Based on the parameters of the remaining two systems and two sets of theoretical 
isochrones obtained for Y=0.25, [Fe/H]=— 1.2 and [a/Fe]=±0.4, we set lower and upper 
limit of the age of M4 at 10.8 and 11.7 Gyr with a formal value of 11.25±0.42 (statistical) 
±0.85 (systematic) Gyr. The isochrone comparison shown in Figure [10] is c onsistent with 
an age ~12 Gyr. An age in excess of 12 Gyr has also been derived for M4 by 



( 120041 ) (from fitting of the white dwarf cooling seque nce; 12.1 Gyr), 



(from CMD fitting based on ACS data; 12.5 Gyr) and iHendricks et al. 



Hansen et al. 



Potter et al. 



(hoioh 



(12012( 1 (from CMD 



fitting based on NTT/SOFI data; 12 Gyr for [Fe/H]=— 1.0 - for a lower metallicity their 
age would be older). The data listed in Table [T41 suggest that this "CMD-DEB discrepancy" 
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mi ght be removed by a slight increase in [Fe/H] sugg ested by the spectroscopic measurements 



of iMarino et al.l (120081 ) and IVillanova et al.l ( 120121 ). It could also be removed by adopting 
modest variations in the fiducial helium content, [a/Fe], or a combination of all three effects. 

We feel, however, that the accuracy of the observational data is still too low, and inherent 
uncertainties in the theory of stellar evolution are still too high, to turn such suggestions into 
firm statements concerning the chemical composition of M4. The first of these two factors is 
illustrated in Figs. [7] and [TU] by the scatter of points which define the main sequence of the 
cluster, and the second - by the sensitivity of stellar evolution codes to det ails of the chemical 
comp osition (see Table HU further discussion of this issue can be found in IVandenBerg et al. 



20121 ) . Formally, considering the statistical and systematic errors, there is no disagreement 
between the two age estimates. 



We note here that lDotter et al.l (120091 ) found a similar result for the DEB V69 in 47 Tuc, 
where the age derived from the properties of the component stars is ~1 Gyr younger than 
the age derived from the location of the main-sequence turnoff (see their Figure 1). Milone 
et al. (2012) have used HST photometry to identify multiple main sequences in the CMD 
of 47 Tuc, concluding that approximately 60% of the cluster population is in the form of 
stars with He and N enrichment. If the ground-based photometry averages out these small 
color differences on the main sequence, then it is difficult to explain the discrepancy in ages 
between the CMD fitting and the age derived from 47 Tuc-V69 as stellar He enrichment since 
the binary appears younger rather than older than the mean pop ulati on, which is already 
appar ently enhanced in He. In the case of M4, IMarino et al.l (120081 ) and I Villanova &: Geisler 
(120111 ) have also found evidence for two populations within the cluster, mainly based on the 
abundances of Na and CN. However there is no structure on the main sequence of the CMD 
that might indicate a He abundance spread. It is not clear how the abundance distributions 
of the two populations might influence age determination through the properties of stellar 
tracks calculated with the different abundances, and we are left with no clear explanation 
of the measured age difference other than assumptions about the chemical composition that 
define the fiducial models used in the age measurement. 

How might the accuracy of the observational data be increased? Given the large incli- 
nations, the errors in the masses of the components of V66 and V69 originate almost entirely 
from the orbital solution, which may only be improved by taking additional spectra (prefer- 
ably with the same instrument). This, however, would require a large observational effort, 
as doubling of the present set of radi al velocity measuremen ts would lead to an improvement 
of only 33% in the mass estimates ([Thompson et al.l 120101 ) . Contributions to the errors in 
the radii of the components are dominated by the photometric solution, whose accuracy, 
in turn, depends on the errors of the differential photometry. The latter originate mainly 
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from a marginally sampled PSF and limited time-resolution, dictated by the diameter of the 
telescopes we used, and the sensitivity of available detectors. We estimate that photometry 
accurate to 0.002 - 0.003 mag in V would reduce the errors of the radii by 50%. Such an 
improvement is entirely viable, as both V66 and V69 reside in relatively sparsely populated 
areas, and are not blended with another stars (see Figure [[]). This goal could be easily 
achieved on a 6-8 m class telescope equipped with a camera capable of good PSF sampling. 
Better data would not remove the anticorrelation of the radii illustrated in Figure |6] and 
briefly discussed in Section E] as this is an inherent property of systems with partial eclipses. 
The axes of the error ellipses, however, would become smaller. On the modeling side, detailed 
evolutionary and atmospheric models made specifically to match M4, for which abundant 
spectroscopic information is available, would improve the accuracy of the age analysis. 

The luminosities of the components are found using absolut e radii and effecti v e tem - 



peratures estimated from (B — V) - T c g calibrations compiled by ICasagrande et al.l ( 120101 ). 
The errors are rather large - in excess of 0.1 x log(L/L ). A significant improvement may be 
expected when IR photometry is obtained, and more accurate calibrations linking (V — K) 
color to surface brightness in V are employed. The anomalous and nonuniform absorption in 
the field of M4 would still have to be accounted for, however the re lation between (V — K ) 



and surface brightness is broadly insensitive to moderate reddening ([Thompson et al" 2001 ). 



Th e agreement of t h e dist ance modulus derived from our DEBS with that recently derived 



by iHendricks et al.l (120121 ) using an entirely different method, together with the fit of the 



observed photometry to model isochrones, suggests that these uncertainties are not too high. 
A further observational test would be to determine the chemical abundances of the 



components of the binaries under study using disentangling software (see e.g. lHadravall2009l ). 
The existing spectra have an adequate S/N to measure velocities but not abundances. Given 
the brightnesses of the components and the orbital periods it is possible to obtain higher 
S/N spectra adequate for abundance analysis purposes. 



We thank the anonymous referee for the detailed and helpful report. This series of 
papers is dedicated to the memory of Bohdan Paczyhski. IBT was supported by NSF grant 
AST-0507325. 
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Table 1. Times of minima for V65 



E HJD-2450000 a 0-C 



1089.0 2402.62228 

1246.5 2763.77804 

1397.0 3108.88034 

1417.0 3154.74157 

1747.5 3912.59451 

2058.0 4624.58352 

2072.0 4656.68639 

2208.5 4969.68786 



0.00024 0.00020 

0.00050 -0.00088 

0.00024 0.00020 

0.00019 -0.00013 

0.00104 -0.00148 

0.00023 0.00018 

0.00027 -0.00005 

0.00093 -0.00079 
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Table 2. Times of minima for V66 



E a 


HJD-2400000 




a 


O-C 


2.0 


49916. 


.64227 


0. 


.00042 


-0.00001 


39.5 


50220. 


,81379 





,00124 


0.00234 


40.0 


50224. 


,87091 


0. 


,00057 


0.00088 


129.5 


50950. 


,83256 





,00172 


0.00089 


174.0 


51311. 


.78618 





.00076 


0.00027 


174.5 


51315. 


.84052 





.00079 


0.00158 


218.5 


51672. 


.73925 





.00032 


0.00020 


219.0 


51676. 


.79534 


0. 


.00045 


-0.00023 


263.0 


52033 


.69286 





.00029 


-0.00040 


263.0 


52033 


,69260 


0. 


.00037 


-0.00014 


308.5 


52402. 


,75682 





,00015 


-0.00006 


308.5 


52402. 


.75675 


0. 


,00025 


0.00001 


353.0 


52763 


,70992 





,00023 


-0.00015 


353.0 


52763 


,71045 


0. 


,00030 


-0.00068 



a Eclipses listed twice were observed in 
both B and V 
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Table 3. Times of minima for V69 



E HJD-2450000 a -C 



87 4240.72872 0.00028 0.00044 
93 4529.85797 0.00032 0.00008 
101. 5 a 4944.69125 0.00096 -0.00042 



a The orbit is eccentric, and the secondary 
minimum occurs at phase 0.6086052(77) 
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Table 4. Equatorial coordinates for three DEBs in M4 (J2000) 



Name 


RA 


Dec 


d a 




h:m:s 


deg:m:s 


arcmin 


V65 


16:23:28.39 


-26:30:22.0 


1.93 


V66 


16:23:32.23 


-26:31:41.3 


0.68 


V69 


16:23:58.01 


-26:37:18.0 


7.69 



a distance from cluster center at RA = 
6:23:35.22, Dec = -26:31:32.7. 
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Table 5. HST photometry for V65 (ID = 13372) and two nearby stars 



Star ID x y V Vega err (V - I) V e 9 a err 



13372 4870.796 4386.209 16.711 0.0022 0.942 0.0031 
13362 4862.341 4393.953 18.728 0.0057 1.053 0.0078 
13374 4877.145 4390.062 17.490 0.0032 0.944 0.0045 
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Table 6. Apparent magnitudes, colors at maximum light, and reddening 



Name 


V 


B — V 


E(B - V) 


V65 


17.028(15) 


0.903(18) 


0.398(10) 


V66 


16.843(12) 


0.878(16) 


0.395(10) 


V69 


17.011(10) 


0.902(10) 


0.403(10) 
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Table 7. Velocity observations of V65 



HJD-2450000 


v p [km s" 1 ] 


v s [km s" 1 ] 


v 3 [km s" 1 ] 


phase 


0780 fiQ778 

z i oz. oy Mo 


~\ ak fin 

140. OU 


3Q sn 
-oy .oU 


R3 39 
00. OZ 


n 0/187 
-U.Z4o ( 


0783 7/in/in 

Z ( OO. ( 4U4U 


k 7n 
-0. 1 u 


1 R7 7/1 
10 I . ( 4 


fiO C\K 
OZ.UO 


n onfin 
u.zuou 


08R8 /1Qfi31 
Zo0o.4y001 


1 7fi 
1.(0 


1 fiO 88 
lOZ.oo 


P.0 KA 
0Z.04 
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62.02 
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4139.85710 


120.45 


5.80 


62.66 


-0.3899 


4259.66316 


129.90 


999.99 


63.01 


-0.1423 


4314.48934 


145.47 


-35.87 


61.96 


-0.2325 


4316.63497 


143.63 


-26.67 


63.99 


-0.2968 


4317.57928 


17.89 


137.38 


62.48 


0.1150 


4317.62312 


11.73 


144.66 


62.49 


0.1341 


4328.48847 


124.59 


-1.92 


62.38 


-0.1275 
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Table 7— Continued 



HJD-2450000 


v p [km s x ] 


v s [km s x ] 


t> 3 [km s x ] 


phase 


4329.49992 


-2.48 


165.31 


62.00 


0.3136 


4966.64890 


-0.07 


162.44 


61.40 


0.1750 


4967.69778 


127.16 


-5.10 


62.14 


-0.3675 


4968.78818 


20.38 


132.90 


61.06 


0.1080 


5012.59393 


-3.46 


171.61 


61.02 


0.2117 


5037.57378 


20.97 


130.66 


61.01 


0.1055 


5354.78703 


37.94 


110.54 


59.04 


0.4426 


5355.65833 


138.62 


-21.08 


61.82 


-0.1775 


5355.70217 


134.42 


-12.90 


62.84 


-0.1583 


5459.51726 


18.60 


140.40 


60.88 


0.1155 
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Table 8. Velocity observations of V66 



TT T T ~\ C\ A rAAAA 

HJD-2450000 


v p [km s J 


v s [km s J 


1 

phase 


2736 81000 


132.56 


20.83 


-0.3163 


2737 79423 


133.76 


18.75 


-0.1950 


2739 75576 


59.89 


95.50 


0.0468 


2867 49150 


136.52 


18.98 


-0 2053 


3066 84849 


36.52 


124.14 


0.3724 


3068 87513 


120.20 


35.70 


-0.3778 


3176.56691 


113.41 


41.07 


-0.1010 


3178.66358 


29.01 


131.01 


0.1575 


3179.53692 


19.63 


141.32 


0.2651 


3179.58097 


19.89 


140.94 


0.2706 


3180.56655 


41.30 


117.88 


0.3921 


3183.53414 


138.30 


16.66 


-0.2421 


3183.57940 


138.06 


16.72 


-0.2365 


3206.63502 


115.50 


40.27 


-0.3941 


3516.81108 


128.63 


28.23 


-0.1541 


3520.65104 


24.98 


135.92 


0.3193 


3584.61769 


22.07 


139.50 


0.2054 


3816.84529 


129.94 


25.27 


-0.1645 


3875.82534 


41.97 


118.00 


0.1069 


3876.68298 


20.77 


139.29 


0.2126 


3877.69367 


27.88 


132.30 


0.3372 
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Table 9. Velocity observations of V69 



HJD-2450000 


w p [km s x ] 


v s [km s x ] 


phase 


3066.80896 


Tr /"* 

75.68 


56.76 


O A O 1 1 

-0.4011 


onpT o o n i r 

3067.82945 


TO O C 

78.25 


ro f 

53.64 


O OTOO 

-0.3799 


3068.82970 


81.34 


r 00 

50.89 


r 00 

-0.3592 


3176.61304 


1 or ^ /i 

105.44 


26.21 


O 1 OO A 

-0.1224 


O 1 T/' 1 f* cr tct c 

3176.65755 


1 Ar A 

l(J5.z4 


oc r n 

2b. 59 


1 01 r 

-U.lzlo 


3178.61436 


00 c 

92.95 


OO CT 1 

39.51 


0000 

-0.0809 


3179.71126 


OO 00 

82.09 


r t r 

50.75 


O O CT O 1 

-0.0581 


Ol on C O /I A f 1 

3182.62446 


53.68 


OO O /I 

80.34 


O OOOO 

0.0023 


oi oo n r\f\ 

3183.62609 


/i T /IT 

47.47 


ot or 1 

87.25 


O OOO 1 

0.0231 


91 no rrr) A fy-\ 

3183.72471 


47.08 


QT T1 

87.71 


n no £ 1 

0.0251 


3184.53725 


AO t K 

43.15 


ai A A 

91.44 


O O A OO 

0.0420 


3201.52003 


r /i nn 

54.00 


TO O A 

79.84 


O OO A A 

U.o944 


oozU. 59751) 


49. Ul 


or 

85.83 


01 OO 

0.0180 


ODzl.OOllZ 


A A 1 /I 

44.14 


nn 

90.75 


n noon 

0.0380 


3581.57763 


/l T1 

44.71 


OO O 1 

90.31 


O OO 1 A 

0.2814 


3582.63902 


An 00 

46.93 


O O O T 

88.87 


O OOO A 

0.3034 


r - /i r-TocrT 

3584.57357 


49.56 


r r\r' 

85.06 


O O A O C 

0.3435 


3585.59952 


CT 1 Ol 

51.81 


O O OO 

83.32 


or* a 

0.3648 


3815.81705 


37.31 


97.99 


0.1423 


3816.75737 


38.34 


97.94 


0.1618 


3817.78611 


38.55 


96.63 


0.1831 


3889.69138 


86.45 


46.25 


-0.3247 


3890.65677 


90.06 


44.13 


-0.3046 


3891.59731 


92.83 


40.37 


-0.2851 


3892.66021 


95.94 


37.23 


-0.2631 


3893.71370 


99.17 


32.86 


-0.2412 


3898.69475 


106.87 


24.08 


-0.1378 


3899.64343 


104.86 


26.98 


-0.1182 


3935.60707 


79.87 


53.09 


-0.3718 


3989.51380 


97.14 


35.16 


-0.2532 
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Table 10. Orbital parameters a 



Parameter 


V65 


V66 


V69 


7 (km s -1 ) 


69.55(15) 


78.76(9) 


66.90(4) 


K p (km s- 1 ) 


77.71(19) 


59.44(14) 


35.28(9) 


K s (km s^ 1 ) 


103.20(50) 


62.62(15) 


37.16(9) 


e 


0.0 b 


0.0 b 


0.3840(12) 


u (deg) 


0.0 b 


0.0 b 


65.25(20) 


(T p (km s -1 ) 


0.99 


0.56 


0.33 


a s (km s _1 ) 


2.54 


0.59 


0.33 


Derived quantities: 








Asmi (R ) 


8.196(25) 


19.561(35) 


63.681(118) 


M p sin 3 % (M ) 


0.8024(86) 


0.7841(45) 


0.7664(44) 


M s sin 3 ? (M ) 


0.6042(44) 


0.7442(42) 


0.7276(43) 



a Numbers in parentheses are the errors of the last signifi- 
cant digit (s) 

b assumed in fit 
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Table 11. Photometric parameters 8. 



Parameter 


V65 


V66 


V69 


i (deg) 


88.30(26) 


89.444(21) 


89.789(11) 


r P 


0.1399(12) 


0.04778(23) 


0.013591(10) 


r s 


0.0745(11) 


0.04242(26) 


0.012681(13) 


(L p / L s ) v 


12.11(28) 


1.489(20) 


1.305(42) 


(L p / L S ) B 


17.54(50) 


1.569(22) 


1.358(44) 


V p (mag) b 


17.114(11)(15) 


17.401(8)(13) 


17.629(15)(18) 


V s (mag) b 


19.822(25)(27) 


17.833(12)(15) 


17.918(20)(22) 


B p (mag) b 


17.991(11)(15) 


18.256(7)(13) 


18.513(14)(17) 


B s (mag) b 


21.101(31)(33) 


18.745(11)(15) 


18.845(20)(22) 


<r ims (y) (mma 


g) 8 


7 


15 


Crms(-B) (mma 


g) 10 


8 


8 


a Numbers in 


parentheses are the 


errors of the last 


significant digits 



b For V p , V s , B p and B s both the internal error (from the photometric 
solution and profile photometry) and the total error is given, the latter 
including 0.01 mag uncertainty of the zero point of the magnitude 
scale. 
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Table 12. Absolute parameters 21 



Parameter 


V65 


V66 


V69 


A (R e ) 




8.200(25) 


19.562(35) 


63.681(118) 


M p (M ) 




0.8035(86) 


0.7842(45) 


0.7665(53) 


M s (M ) 




0.6050(44) 


0.7443(42) 


0.7278(48) 


R p (R Q ) 




1.1470(104) 


0.9347(48) 


0.8655(97) 


i? s (R Q ) 




0.6110(92) 


0.8298(53) 


0.8074(80) 


r P (K) 




6088(108) 


6162(98) 


6084(121) 


r.(K) 




4812(125) 


5938(105) 


5915(137) 


^ bo1 (L ) 




1.620(118) 


1.129(73) 


0.920(76) 


(Lo) 




0.179(19) 


0.767(55) 


0.715(68) 


log[g- p (cm s 


- 2 )] 


4.221(14) 


4.388(77) 


4.444(132) 


log[g- s (cm s' 


- 2 )] 


4.645(17) 


4.469(85) 


4.483(119) 


M v P (mag) 




4.329(82) 


4.716(74) 


4.945(91) 


My s (mag) 




6.983(114) 


5.153(81) 


5.230(103) 


(m - M) Vp ( 


mag) 


11.400(94) 


11.310(88) 


11.281(103) 


(m - M) Vs ( 


mag) 


11.454(110) 


11.305(90) 


11.285(112) 



a Numbers in parentheses are the errors of the last signifi- 
cant digits 



Table 13. Results from Isochrone Age Analyses 



ID Mass-Radius Mass-Luminosity 
Age (Gyr) Age (Gyr) 



Dartmouth 



V65A 


12.087 ±0.712 


11.618 ±0.998 


V65B 


15.068 ± 1.968 


10.126 ±4.033 


V66A 


10.763 ± 0.442 


10.722 ±0.861 


V66B 


11.496 ±0.587 


11.217 ± 1.288 


V69A 


10.570 ±0.613 


10.829 ± 1.277 


V69B 


12.382 ±0.779 


12.677 ± 1.701 




Victoria-Rej 


$ina 


V65A 


11.797 ±0.763 


11.639 ± 1.000 


V65B 


16.330 ± 1.035 


10.428 ±2.104 


V66A 


10.617 ±0.491 


10.948 ± 0.835 


V66B 


11.638 ±0.605 


11.769 ± 1.215 


V69A 


10.587 ±0.657 


11.228 ± 1.215 


V69B 


12.595 ±0.800 


13.224 ± 1.613 



Note. These results assume 

[Fe/H] = -1.2, [a/Fe]=+0.4, and Y=0.25. 
All ages are given as mean ± standard 
deviation derived from the age histograms 
presented in Figure [9] (see text for discus- 
sion) . 
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Table 14. Age Sensitivity of the Mass-Radius Plane to Chemical Composition 



[Fe/H] 


[a/Fe] 


Y 


A- age 


Model 








(Gyr) 




-1.0 


+0.4 


0.25 


0.9 


Dartmouth 


-1.2 


0.0 


0.25 


-0.4 


Dartmouth 


-1.2 


0.0 


0.25 


-1.0 


Victoria-Regina 


-1.2 


+0.4 


0.27 


-1.7 


Victoria-Regina 


-1.2 


+0.4 


0.29 


-3.2 


Victoria-Regina 



Note. — A-age is calculated such that a positive 
value means that the model with varied composi- 
tion yields an older age than the model with the 
fiducial composition ([Fe/H] = — 1.2, [a/Fe]~ +0.4, 
and Y=0.25). 




Fig. 1. — Finding charts for V65 (left), V66 (middle) and V69 (right). Each chart is 26" on 
a side, and is oriented with north up and east to the left. The coordinates of the targets are 
given in Table HI 
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Fig. 2. — Complete photometric observations of V65, illustrating the instability of the light 
curve. The median internal errors in B and V are, respectively, 0.009 mag and 0.007 mag. 
Note that to bring the data closer together on the plot the 5-curve has been shifted to lower 
values by 0.7 mag. 
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Fig. 3. — Light curves adopted for the analysis: V65 (left), V66 (middle) and V69 (right). 
The light curve of V65 only includes observations collected between 2008 June 07 and 2009 
June 30. Note that to bring the data closer together on the plot the 5-curves have been 
shifted to lower values by 0.8 mag (V65) and 0.7 mag (V66, V69). 
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Fig. 4. — Velocity curves adopted for the analysis: V65 (left), V66 (middle) and V69 (right). 
Since the spectrum of V65 was contaminated by one of the stars blended with the target, the 
measurements were performed with a three-dimensional extension of the original TODCOR 
algorithm. The velocities of that star are marked in the left panel with triangles. 
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Fig. 5. — Residuals to the fits to the light curves adopted for the analysis: V65 (top), V66 
(middle) and V69 (bottom). In each panel, the lower sequence represents the l^-residuals, 
and the upper one - the 5-residuals offset by 0.1 mag for clarity. 
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Fig. 6. — Relative radii of the components in V65 (left), V66 (middle) and V69 (right). White 
crosses: light curve solutions from Table [TTJ Points: results of Monte Carlo simulations 
described in Section HJ Inner contour: confidence level 67%. Outer contour: confidence level 
93%. Whenever the errors were not symmetric, the larger one was listed in Table [Til an d 
used for further error analysis. 
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Fig. 7. — The observed color-magnitude diagram of M4 with locations of the components of 
the three systems investigated in the present paper. Data are taken from Kaluzny et al. (in 
preparation) . 
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Fig. 8. — Radius-mass and luminosity-mass diagrams with locations of the components of 
the three systems investigated in the present paper. Solid and dotted lines are Dartmouth 
and Victoria-Regina isochrones, respectively, for 13, 12, 11 and 10 Gyr from left to right. 
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Fig. 9. — Age histograms based on the mass-radius (top) and mass-luminosity (bot- 
tom) analysis for Dartmouth and Victoria-Regina isochrones, both assuming [Fe/H] = — 1.2, 
[a/Fe]~ +0.4, and Y=0.25. Heavy lines: V66 primary (solid) and V66 secondary (dotted). 
Thin lines: V69 primary (solid) and V69 secondary (dotted). 
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Fig. 10. — The observed CMD of M4 (the same as plotted in Figure [7j) compared with 
Dartmouth isochrones calculated for the fiducial composition [Fe/H] =— 1.2, [a/Fe]~ +0.4, 
Y=0.25, and ages 10, 11, 12, and 13 Gyr (the same as plotted in the mass-radius and mass- 
luminosity planes in Figure [8]). PHOENIX synthetic fluxes were used, and the isochrones 
were adjusted for a true distance modulus 11.34, E(£> — V)=0.39 and v4y=1.47. 



